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Abstract 

Lorentz-invariant scalar field theories in d + 1 dimensions with second-order derivative terms are 
unable to support static soliton solutions that are both finite in energy and stable for d > 2, a result 
known as Derrick's theorem. Lifshitz theories, which introduce higher-order spatial derivatives, 
need not obey Derrick's theorem. We construct stable, finite-energy, static soliton solutions in 
Lifshitz scalar field theories in 3 + 1 dimensions with dynamical critical exponent z = 2. We 
exhibit three generic types: non-topological point defects, topological point defects, and topological 
strings. We focus mainly on Lifshitz theories that are defined through a superpotential and admit 
BPS solutions. These kinds of theories are the bosonic sectors of super symmetric theories derived 
from the stochastic dynamics of a scalar field theory in one higher dimension. If nature obeys a 
Lifshitz field theory in the ultraviolet, then the novel topological defects discussed here may exist 
as relics from the early universe. Their discovery would prove that standard field theory breaks 
down at short distance scales. 
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I. INTRODUCTION 



Hofava has suggested that nature may exhibit Lifshitz-like anisotropic scaling in the 
extreme ultraviolet (UV) regime in order to cure the divergence problem in quantum gravity 
l|, l2j . If gravity were to behave this way, then presumably all other interactions should also. 
In this paper, we shall show that Lifshitz scalar field theories can support stable, finite- 
energy soliton or defect solutions of a kind that are impossible in standard second-order 
Lorentz-invariant field theories. Specifically, we shall demonstrate that stable, static, finite- 
energy defects formed entirely from scalar fields can exist in three dimensions. We shall 
derive both point-like and string-like solitons, which we can call "Lifshitz poles (L-poles)" 
and "Lifshitz strings (L-strings)" , respectively. The non-existence of solutions having all 
of the above characteristics in standard field theories follows from Derrick's theorem 
So, while point-like ("hedgehog") and global string solutions do exist in standard theory, 
they have divergent energies or linear energy densities, as the case may be. The discovery 
of L-poles and/or L-strings as cosmological relics from the early universe would provide 
observational proof of the breakdown of standard field theory at short distances. 

The basic idea of Lifshitz theories is to add high-order spatial derivative terms to the ac- 
tion while maintaining time derivatives at second order. This procedure preserves unitarity, 
while softening UV divergences through the presence of higher powers of 3-momentum in 
the denominators of propagators. This is an interesting idea, because it increases the range 
of power-counting renormalisable field theories. Hofava gravity jl| is the best-studied ex- 
ample; field-theory models involving extra dimensions of space is another [4|. The improved 
UV behaviour comes at the expense of Lorentz invariance, with the UV isotropic scaling 
of Lorentz-invariant field theories replaced by anisotropic scaling between space and time. 
This is obviously the most serious drawback of this approach, in both a phenomenological 
and an aesthetic sense. It is not known how to ensure that standard relativity emerges in 
the infrared, except by fine-tuning, though progress has recently been made on the gravity 
front 2]. More generally, it is neither clear how to ensure that all particle species acquire a 
common limiting speed, nor how to make the renormalisation group running of the limiting 
speed(s) phenomenologically acceptable J5]. So, there are arguments both for and against 
Lifshitz theories of particle physics and gravity. Experimental or observational evidence in 
favour of the breakdown of Lorentz invariance in the UV would lend support to the Lifshitz 
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approach. Such evidence can be searched for in very high energy cosmic ray events, but also 
through the possible existence of cosmological relics that would be inexplicable by standard 
field theory. 

We shall examine Lifshitz theories having fourth-order spatial derivatives, which is the 
simplest possible choice from the point of view of the differential equations. In general, 
these theories are much more difficult to analyse than the usual second-order relativistic 
theories. However, a Lifshitz analogue of the Bogomolnyi-Prasad-Sommerfield (BPS) 6] 
superpotential approach can be invoked to reduce the problem back to second-order (this 
is equivalent to the "detailed balance" condition |7J). We shall take this approach in this 
paper, partly for simplicity, and partly because the BPS system has intrinsic theoretical 
interest through a deep connection with supersymmetry and stochastic quantisation. 

In the next section we review Lifshitz scalar field theory and anisotropic scaling. We then 
show why Derrick's theorem is inapplicable. This is followed by the derivation of Lifshitz 
BPS soliton solutions, and then a brief section on the general fourth-order system. The final 
section is a conclusion. 



II. LIFSHITZ SCALAR FIELD THEORY 



We shall work exclusively in 3 + 1 dimensions. The construction of a Lifshitz theory 
begins through the specification of the dynamical critical exponent z, 

X -> KX , t -> K Z t, (2.1) 

which defines the anisotropic scaling invariance of the theory in the UV. 1 Lorentz- invariant 
theories must have isotropic scaling, namely 2 = 1. Consider now the case of z = 2, adopted 
for definiteness and simplicity in our analysis. We shall construct the action step-by-step. 
One first chooses that the Lagrangian contains only the term <p 2 in time derivatives. This 
is to ensure the unitarity of the quantum theory. Anisotropic z = 2 scale invariance then 
follows with the observation that the action J d 3 x dt </> 2 is invariant under 

x — y nx, t -> kH, 4>{x,t) -» n~ 1/2 <f)(x,t), (2.2) 

1 We shall not consider the possibility that different dimensions of space may scale differently. 
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or, equivalently, under the field substitution 

t) -> k +1/2 <P{kx , kH), (2.3) 

without the coordinates in the integration measure and derivatives participating in the 
transformation. We shall adopt the viewpoint of Eq. (I2.2p for definiteness in this paper. 

The 2 action is actually invariant for any z, with the field transforming as — > n^cj). 
The z value is specified through the spatial-derivative terms we now introduce. The rela- 
tivistic choice, J d 3 xdt (4> 2 — c 2 V0 • V0)/2, enforces z = 1. But we are interested in z = 2, 
which implies that 

S C J d z xdt l - [</> 2 - (V 2 0) 2 ] . (2.4) 

The arbitrary constant in front of the second term has been absorbed into x, and the minus 
sign is necessary to make the energy bounded from below. 2 The higher spatial derivative 
implies that the propagator behaves like l/[k 2 — (k ■ k) 2 ], so has improved UV behaviour 
compared to the relativistic case 1/(&q — k ■ k). 

The most general z = 2 scale invariant action for a single real scalar field is then easily 
deduced to be 3 



S = d 3 x dt 



^ 2 -^(VW-a0 5 V 2 0-A0 10 



(2.5) 



where a and A are coupling constants, and they are dimensionless with respect to the 
"anisotropic scaling units", 

[S] = 0, [x] = -l, [t] = -2, [0] = i (2.6) 

that replace the usual "mass dimension units" of relativistic theories. The power-counting 
renormalisation properties follow from the adoption of these new units, with the usual rule 
that renormalisable terms have either dimensionless or positive-power coupling constants, 
corresponding to marginal and relevant operators, respectively. The scale-invariant action 
( 12. 5ft contains only marginal operators. The most general power- counting renormalisable 
action contains in addition the terms, 

n V 2 0, with < n < 4 (2.7) 
(p m , with < m < 9. (2.8) 



2 Note that the constant must be reintroduced if one wishes to work in standard natural units. 

3 Note that we are assuming that surface terms may be discarded, so that there is no need to include 
4 V0 ■ V<jf>, V0 ■ V(V 2 </>) and <^>V 2 (V 2 0) as independent terms. 
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These relevant operators explicitly, but softly, break the anisotropic scale invariance. The 
most general action may be compactly written as 



S 



d xdt 



if(0)V 2 - V((f>) 



(2.9) 



where K(4>) is quintic and the potential V(4>) is tenth-order. Observe that the standard 
relativistic term V0 • V0 is contained as the n = 1 term of K, up to a discarded total 
derivative. The equation of motion that follows from Hamilton's principle is 



+ V 2 (W) - K H V(j) ■ V0 - 2AT VV + V+ = 0, 



(2.10) 



where the subscripts and 00 denote first and second derivatives with respect to 0, respec- 
tively. The energy functional is 



E[4>] 



d 3 x 



I02 + I(v 2 0) 2 + ^(0)V 2 + V(<f>) 



(2.11) 



The generalisation of the above formalism to more than one scalar field is obvious. 



III. DERRICK'S THEOREM AND THE VIRIAL RELATION 



Derrick proved that stable, finite-energy, static solutions to second-order relativistic scalar 
field theories do not exist in greater than two spatial dimensions We now review his 
extremely simple argument for the case of three spatial dimensions. For static configurations, 
extremising the action functional is equivalent to extremising the energy functional, 



E[<f>] = J d 3 x 
Consider the specific field variation given by 



-V0-V0 + V(0) 



under which the energy varies to 

£,.M = / d 3 x 



4>(x) — > 4>(kx), 



-V(f){kx) ■ V(j){kx) + V {(j>{kx)) 



(3.1) 



(3.2) 



(3.3) 



A necessary condition for a solution is that the variation of the energy under this restricted 
field variation should vanish to first order, which is equivalent to the demand that 

dE k k 



dk 



0. 



(3.4) 



k=l 



Now 



E h [<f>\ = 



j d?i fc-^v f 0(e)-v f 0(o+^v(0(e)) 

AT 1 Si + k~ 3 S 3 



(3.5) 



where £ = kx and E — Si + S 3 . Imposing Eq. (13. 4p we get the virial relation — Si = 3S 3 , 
which implies that E = 2S\/3. A stable solution must have a positive second-order variation. 
Computing d 2 E k /dk 2 \ k=1 one obtains 2 Si + 12S 3 , which equals — 2 Si through the virial 
theorem. Since Si is positive-definite, no stable solutions can exist. It was implicit above 
that E[(f)] converged, so the theorem does not preclude infinite-energy, static, point-like 
solitons (the hedgehog being a well-known example). 

Repeating this exercise for the Lifshitz energy functional (12.111) we obtain 



where S± = J d 3 xdt (V 2 0) 2 /2 and S2 = J d 3 x dt K((p)'V 2 <j), which implies the virial relation 



The second-order variation is now equal to 2S 2 + 12S 3 , which in turn equals 4Si — 2S 2 
through the virial relation. There is no necessity for this quantity to be negative, hence 
there is no generalisation of Derrick's theorem to the z = 2 Lifshitz case. A simple example 
is the case K — which implies S 2 = and hence a positive second-order variation. Another 
is K = (p 5 , which sees S 2 turn into an integral of the negative-definite function — 4 V0 • V0 
after integration by parts. This makes 4Si — 2S 2 necessarily positive. The addition of higher- 
derivative terms thus tends to stabilise scalar-field solitons, a phenomenon also observed in 
the (relativistic) Skyrme model [9(. However, evasion of Derrick's theorem does not prove 
that stable, static, finite-energy, scalar- field solitons actually exist, because the variation 
(13. 2p is not general. We therefore turn to explicit computations. 

IV. THE LIFSHITZ BPS SOLITON SOLUTIONS 

A technical problem with z = 2 Lifshitz theory is that the scalar field equations are fourth- 
order in spatial derivatives, as per Eq. (I2.10p . To get to concrete results in the simplest 
possible way, we therefore develop a Lifshitz generalisation of superpotential systems. This 



E k [(j)\ = kSi + A; _1 S 2 + k' 3 S 3 , 



(3.6) 



Si — S 2 — 3S3 — 0. 



(3.7) 
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will reduce the problem to second order, and as a bonus the solutions will be Lifshitz 
generalisations of BPS states, which is interesting in itself because of the strong connection 
with supersymmetry and stochastic quantisation. We shall derive three kinds of solutions: 
a perturbatively-stable non-topological point-like soliton, a topological point-like soliton or 
hedgehog, and a topological string. In the section after this, we shall also briefly discuss the 
general fourth-order problem, and display a finite-energy solution by way of an existence 
proof, but defer a systematic treatment to future work. 

We first very briefly review how superpotential considerations lead to interesting simpli- 
fications in one-dimensional second-order systems, with 



E\ 



dx 



For the special cases where 

with W being a superpotential, it is easy to see that 



E\ 



=F [W {<j>{x = +00)) - W (<f>{x = -00))] + 



dx ( — ± 

dx 



dW 



Since the integrand above is positive definite, solutions to the first-order equations 

dW 



dx 



± 







(4.1) 



(4.2) 



(4.3) 



(4.4) 



globally minimise the energy functional for given boundary conditions. They are stable BPS 
solutions, and it is trivial to confirm that they satisfy the Euler-Lagrange equations. 

We can achieve a situation that resembles the above if we restrict our Lifshitz action 
functional to the special form 

,2 



S 



d 6 xdt 



so that 



1 



v{4>) = ^(0) 2 . 



(4.5) 



(4.6) 



A superpotential W can then be defined as the negative indefinite integral of K with respect 
to <ft, so that K = —Wff) (the minus sign is just a convention). The convenience of the action 
(14. 5 p is that the energy functional, 

, 2 



d x 



\^ + \{v 2 <p + K(<p)y 



(4.7) 



is positive definite. Configurations that obey 

= 0, V 2 (j) + K((f)) = 



(4.8) 



minimise the energy, and in fact they have zero energy- density. It is straightforward to 
confirm that solutions of Eq. (14.81) are also solutions of the field equation (I2.10p . We shall 
call these "Lifshitz BPS" or LBPS solutions. 

From Eq. (12.101) . we see that static, spatially homogeneous solutions (vacua <fi = v) 
must obey V^{v) = as usual. For a Lifshitz superpotential system, this means that 
K{v)K(f,(v) = 0, which implies that either K(v) = 0, or K^{v) = with K(v) ^ 0. The 
former type of vacuum obeys Eq. (I4.8p . while the latter does not. Suppose K{<p) ^ V0. 
Then LBPS solutions must obey either V 2 > everywhere, or V 2 < everywhere. 
In either case, these solutions will be asymptotically unbounded, and thus forbidden once 
physically reasonable boundary conditions are imposed. We shall not consider such theories 
any further, by requiring that K vanishes for all vacuum configurations. This means the 
LBPS solutions for the allowed models are degenerate with the vacua, and hence at least 
perturbatively stable. 

An interesting aspect of Lifshitz theories with z = 2 anisotropic scaling (14.51) is that they 
provide an equivalent description of certain stochastic systems Q, [ll| (for recent relevant 



works see [12|] and references therein) . It is well known that the equilibrium limit — defined 
as the r — > oo limit where r is a fictitious time (Markov parameter) — of the statistical 
correlators of a (d+ l)-dimensional scalar field 0^(r, x l ) is equivalent to quantum correlators 
of the d-dimensional 4 Euclidean scalar field (f)(x l ) with the classical action S*e [</>]: 

lim (0„(r, x\)..4r,{r, <))„ = (<j>(x{)..4(j n )). (4.9) 

T— J-OO 

Here 4> v (t, x l ) is a stochastic field, being a solution of the Langevin equation, 



a r <f>( T ,X i ) = -^+Tl(T,X i ), (4.10) 



where 



SSk 5Sk 



(4.11) 

(x 1 )=4>(t,x' 1 ) 



In what follows we explicitly discuss the d=3 case. 
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and rj(r,x l ) is Gaussian white noise, i.e., 



(t]{t', x h )t]{t, x l )) v = 5{t' - t)5 3 {x h - x l ) 



(4.12) 



with all other correlators being zero. Note that in the large time limit, r — > oo, the system 
approaches equilibrium with <9 t 0(t, x l ) = r\ = 0, and the Langevin equation (I4.10p . reduces 
to the equation of motion of the Euclidean scalar field <fi(x l ) described by the action Se- 
The generating functional for the statistical correlators in Eq. (14. 9 p reads as: 

,2 



Z[J]= / Drjexp { / drd 3 x 



and 



1] 

■~2 + J 0v 



5 n Z 



(4.13) 



(4.14) 



j=o 



Z 6J(T,ai)...6J{T,x\) 

Following the standard steps, we perform change of variables r](T,x l ) — > 4> v { T i x l ) using 
Eq. (I4.12p and express the new functional measure through a path integral over the Grass- 
mannian fields ip and ip: 



Dr/ ->• D(p det 



D<j> det d T + 



5 2 S B 



D<f> / DipDipexp < / drd 3 xtp d T + 



5 2 Sv 



(4.15) 



Plugging (I4.15P into (I4.13P and using the Langevin equation, we obtain the effective action 

^ 2+ Hf) 2 -^ + 5 N ' • (416) 

where we have dropped the full time derivative term <9 T 0^rf . Taking 



S'cff = / drd x 



ft 



d 3 x 



l -(d % <P) 2 + W{<P) 



(4.17) 



in (I4.16P and passing to the Lorentzian time t = —it, we obtain a theory described by the 
action: 

2 / 721 



s 



LSUSY 



dtd 3 x 



' ^- • 4 + i/j ( -id t - V 2 + ) V 



(4.18) 



The bosonic part of the above action is precisely the action of the z = 2 Lifshitz theory (14. 5 j) 
with = -W . 



The action (I4.18P is invariant under rigid (quantum mechanical) supersymmetry trans- 
formations. To obtain transformations that are linear in the fields, we rewrite f 14 . 1 8 1) in the 
equivalent form 



Slsusy = / dtdTx 



where D(t,x l ) is an auxiliary field satisfying the non-dynamical equation of motion: 

D = V 2 (j) _dW_ (420) 
d<p 

The infinitesimal supersymmetry transformations then read: 

5 e <fi = eip - $e, 5 e D = -id t (ei/j + ipe), 
8 e ip = e(D - ij>), 5^ = -e(D + #). (4.21) 

It is clear now that any solution of z = 2 Lifshitz theory that satisfies Eq. (14.81) gives 
vanishing D-term (I4.20p . Therefore, our LBPS solutions are supersymmetric and, hence, 
they are degenerate with the zero-energy ground state. Supersymmetry ensures the stability 
of the LBPS solitons. 

From the perspective of stochastic dynamics (14.101) . our supersymmetric LBPS solitons 
are relevant to the description of physics at equilibrium (r — > oo), while any time- dependent 
solution describes physics in the non-equilibrium regime. A discussion of these interesting 
points is beyond the scope of the present paper. 

A. Non-topological BPS L-pole 

Our first explicit solution will be a non-topological BPS L-pole. Let = <po(r) be a 
radially symmetric solution of Eq. (14. 8ft : 

+ (.Mr))- (4-22) 

Impose the boundary condition that (fro asymptotes to a finite value v as r — > oo. Regularity 
at the origin requires 0o ~ r2 as r — * 0- An obvious interpolation between these behaviours 
is provided by the function 

Mr) = (4-23) 
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FIG. 1: Numerical solution for the non-top ological BPS L-pole corresponding to the choice of K 
given in Eq. f)4.25j) with A = v = 1. 

where R is a parameter that sets the characteristic size of the defect. It is easy to verify 
that this function solves Eq. (14.221) for 

R z v z 

so that v is a vacuum. This provides a simple analytical example of an L-pole. It is non- 
topological, because the asymptotic limit is always the same vacuum. 

To get the simple analytical form of Eq. (I4.23P we of course need a very specific function K. 
If we vary this function continuously from Eq. (I4.24p . then the solution will vary continuously 
away from our analytical example, with numerical computations being in general necessary. 

We give another example of a non-topological BPS L-pole which corresponds to a simpler 
choice of K, namely 

K(<f>) = A0(0 2 - v 2 ) 2 . (4.25) 

The non-topological BPS L-pole corresponding to this choice of K does not have an analytic 
solution, however a numerical solution is shown in Fig. [TJ 
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r 



FIG. 2: Numerical solution to Eq. (g3gj) with A = v = 1. 



B. Topological BPS L-pole or hedgehog 



Our second example is an LBPS hedgehog. The goal is to produce a topologically non- 
trivial BPS L-pole. We need the vacuum manifold to be a 2-sphere, so we choose an 0(3) 
invariant theory with a triplet of scalar fields = (0j), % = 1,2,3, and work within the 
parameter regime that features 0(3) — » 0(2) spontaneous symmetry breaking. The energy 
functional for static configurations is 



E[<f}= J d 3 x 
where there is a sum over i and 



V 2 ^ + <p t F{<p ■ 4>) 



F = A(0-0-w 2 ) 2 . 
The hedgehog ansatz in spherical polar coordinates is 

<f>i(r, 9,<p) = Xih(r) 



(4.26) 



(4.27) 



(4.28) 



where (xi, x 2 , x 3 ) is a unit vector in the (9, <p) direction. The LBPS equations V 2 0j+0jF = 
reduce to 



d 2 h 2dh h 



+ 



dr 2 r dr r 



2- + Xh(h 2 -v 



2\2 



0. 



(4.29) 
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which is the same as the original hedgehog equation jg] except that the non-derivative term 
is fifth-order rather than third-order. We demand that it asymptote to v and be regular at 
the origin. A numerical solution is depicted in Fig. [2J 



C. Topological BPS L-string 

The third solution is a topological BPS L-string. The theory has a complex scalar field 
$ and a spontaneously broken $ — > e ia Q global U(l) symmetry, with energy functional: 

E[$] = J d 3 x [V 2 <F + A$*($*$ - v 2 ) 2 ] [V 2 $ + A$($*$ - v 2 ) 2 ] . (4.30) 

The topological global string ansatz is 

*(p, 0) = e"/(p) (4.31) 

where (p, 6, z) are the usual cylindrical coordinates, and the string defines the z-axis. The 
BPS equation of motion V 2 $ + A$($*<1> - i; 2 ) 2 = reduces to 

d 2 f ldf f 



(4.32) 



dp 2 p dp p 2 

The asymptotic boundary condition is / — > v as p — > oo, and we demand regularity at the 
origin. A numerical solution to Eq. ( 14.32}) is given in Fig. |3j 




FIG. 3: Numerical solution to Eq. (TQ2]) with \ = v = l. 
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V. GENERAL LIFSHITZ SYSTEMS 



One may wonder if the existence of the finite-energy, static solutions displayed above relies 
on the special superpotential form of the action. The answer is no, which we shall establish 
by displaying a numerical solution for a non-LBPS hedgehog. A full analysis of non-LBPS 
solutions is beyond the scope of this paper, in part because the stability of non-topological 
soliton solutions requires the development of a fourth-order linear stability formalism. 

We shall exhibit a non-LBPS hedgehog solution for an 0(3) triplet model with K = 
and 

V=±($-$-v 2 ) 5 , (5.1) 
where, using the ansatz of Eq. (14.281) . the equation of motion reduces to 

^ + 7*1-^ + ^-^ = 0. (5.2) 

A numerical solution for h(r) is given in Fig. fll]). This configuration is topologically stable. 




01 23456789 10 



FIG. 4: Numerical solution to Eq. (I5.2p with X = v = 1. 
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VI. CONCLUSIONS 



Lifshitz scalar field theories permit the existence of finite-energy, static soliton solutions 
that have only infinite-energy analogues in standard second-order relativistic theories (such 
as the global string and the hedgehog). We have displayed several such solutions, mostly 
within a superpotential-based subclass of z = 2 Lifshitz scalar theories that admit stable 
zero energy density BPS solutions. These systems constitute the bosonic sector of super- 
symmetric theories that are themselves related to a stochastic field theory in one higher 
dimension. If the world is described by a Lifshitz theory, as suggested recently by Hofava 
as a way to improve ultraviolet behaviour in a ghost-free manner, then the existence of 
new classes of finite-energy defects may provide a way to look for observational evidence in 
support of this hypothesis. To really make this connection, one would need to extend our 
analysis from z = 2 (which was chosen for simplicity), to at least z = 3. The gravitational 
properties of Lifshitz solitons could then also be studied. 
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